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Abstract 

The spectrum of excitations of the chiral superconducting ring with internal and external radii 
Ri, Re (comparable with coherence length ^) trapping a unit flux $0 is calculated . We find 
within the Bogoliubov-deGennes approach that there exists a pair of precisely zero energy states 
when 2k±_ [Re — Ri) /vr is integer (here k± is the momentum component in the disk plane while 
k±^ > 1 ). They are not protected by topology, but are stable under certain deformations of the 
system. We discuss the ways to tune the system so that it grows into such a "Majorana disk". 
This condition has a character of a resonance phenomenon. 
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I. INTRODUCTION. 



Spin-triplet p-wave superfluids, both neutral, such as liquid He^, Li^ or and charged, 
such as superconducting material Sr^RuO^ and possibly heavy fermion UPt^ have resulted 

n 

in very rich physics [1]. The condensate is described by a generally tensorial complex order 
parameter A exhibiting a great variety of the broken symmetries ground states. The broken 
symmetry and boundary conditions give rise to the continuous configuration of the order 
parameter as nontrivial topological excitations {2!. Especially interesting is the case of the 
so-called topological superconductors, characterized by the presence of electron-hole sym- 
metry and the absence of both time-reversal and spin-rotation symmetry. Realizations of 
topological p-wave superfluids are chiral superconductors like Sr2RuOi, with order param- 
eter of the Px ± ipy symmetry type sl and the ABM - phase jl| of superfluid He^ and other 
fermionic cold atoms, as well as the topological superconductor CuxBi2Se3 that produces 
an equivalent pseudospin system on its surface j4|. 

Generally a magnetic field in type II su|)erconducting films easily creates stable line - 
like topological defects, Abrikosov vortices[5!]. In the simplest vortex the phase of the order 
parameter rotates by 27r around the vortex and each vortex carries a unit of magnetic 
fiux $0 with superfiuid density depleted in the core of the size of the coherence length ^. 
Quasiparticles near the vortex core "feel" the phase wind by creating a set of discrete low- 
energy Andreev bound state. An unpinned vortex in an infinite p-wave superconductors 
exhibits a remarkable topological feature of appearance of the zero energy mode in the 
vortex core jsl (for each value of momentum k± perpendicular to the field direction). The 
zero mode represents a condensed matter analog of the Majorana fermion j^]. Its topological 
nature ensures robustness against perturbations from deformations of order parameters and 
nonmagnetic impurities. Due to possible applications of the Majorana states in quantum 
computing it . important to en.u.e a .datively large min.gap E^, separating the Majorana 
states from charged excitations. It was proposed to enlarge the minigap from [8|] A'^/Ep to 



lo| i?^ ~ e 



order A by pinning vortices on inclusions of small radius 

In this note we study the infiuence of size effects on the appearance of the Majorana 
states in a ring made of the chiral superconductor. It was noted that in finite samples of size 
L » ^ there generally is present there is a pair of "nearly" Majorana states constituting 
one charged fermion degree of freedom like in ID Kitaev model ll[ [l2|. However these states 
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are no longer topologically protected and naively one would expect that the degeneracy is 
removed with small splitting energy 5E cc e~^^^ 13|, Uj^ due to tunneling from the core to 
the surface of the sample. For systems of small size L ~ ^ the situation might be different. 
In particular it is not clear whether the splitting always occurs at all and how it depends of 
the sample geometry. 

To investigate the survival of Majorana states we calculate exactly the spectrum of ex- 
citations in the chiral superconducting ring with internal and external radii Ri, Re of order 
^ trapping a unit flux $o (see fig.l) within the Bogoliubov - deGennes approach. A sur- 
prising finding is that although the splitting is generally of order of the bulk energy gap A, 
there exists a pair of precisely zero energy states for special relation R^, — Ri ^ Tin/2k_^_ for 
any integer n. Then, using perturbation theory, we generalize the consideration to other 
geometries. 



II. BASIC EQUATIONS. 

We start with the Bogoliubov - deGennes (BdG) equations for the Px + ipy superconductor 
in the presence of a single pinned vortex. The vector potential A in polar coordinates r, ip 
has only an azimuthal component (r) and in the London gauge consists of the singular 
part = hc/2er and the regular part of the vector potential that can be neglected for a 
type II superconductor jl5|. In the operator matrix form for a two component amplitude 
the BdG equations read: 



Ho L 

L+ -H* 





(1) 



where for anisotropic dispersion 



Ho = 


vi 




L = 





2m, 



- Er, (2) 



with A being the "bulk gap" of order Tc (neglecting the small inhomogeneity of the superfiuid 
density within the ring). The dimensionless profile of the order parameter s (r) is defined 
to represent the gap function A (r) = As(r). In principle it should be determined self 
consistently, however, for a sufficiently thin homogeneous disk we initially take s (r) = 1. 



This is justified in the bulk since the sample size is of order of ^. Although the Andreev's 
bound states are typically inhomogeneous, the effect of their inhomogeneity on the order 



parameter is still smaller than that of the continuum states even for small sizes ^\. We will 
later investigate the stability of the solutions with respect to variations to s (r). 

The equations possess the rotational and the electron-hole symmetries and eigenstates 
can be found in a form: 



u = (r) e'^^e'''"' 
V = (r) e*('-2)V="^ 



(3) 



For any angular momentum / and momentum along the field k„, there are radial excitation 
levels. In a dimensionless form Eq.([T]) is written as 

_ + i A _ i! + f _ (2— - g - e f- (4) 

^ /p2 47^ / y Ot T J 

with the dimensionless energy eik„ = EikJ'yA. Here distances are in units of ^. In the clean 
limit BCS (applicable to SrRu204) C, = hk^/m^A, where 

kl/2m± = Ep/tf - A;2/2m,„ (5) 

and for given k„ there is just one dimensionless parameter 

7 = l/2k^^ = m^A/2h'^kl « 1. (6) 

The Ansatz, Eq.Q was chosen in such a way that the equations become real. In a 



microscopic theory of the superconductor-insulator interface, (see [16|), the order parameter 
rises abruptly from zero in a dielectric, where amplitudes of normal excitations / (Re) = 
g {Ri) = 0, to a finite value inside the superconductor within an atomic distance a from 
the interface, namely with a slope oc 1/a. This means that the boundary condition on 
the amplitudes is consistent with a zero order parameter at the boundary point in the 
self-consistency equation (see details in ref. jo] and references therein). 

III. MAJORANA VERSUS NON-MAJORANA RINGS. 

Let us first determine under what conditions zero energy (Majorana) states appear. It 
was shown 1^ that they appear only for / = 1 and might contain two possible states 



(7+ = /+ or/and = — /_. The corresponding equations simplify and differ just by the 
transformation r — — r: 

— + (--2] — ---- + — ^ / = 



The solutions are 



f± = e 



c±Ji ( -rJ-^-l) +d^Yi (rJ-^-1 



(8) 
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where Ji and are the Bessel functions. Defining the ratio x (^) = 

Ji ^ry^^^5— /Fi ^— ry^^ — plotted in fig. 2, the boundary conditions for both /+ 
and /_ read: 

X (Ri) = X (Re) = -d±/c± = X (9) 

Therefore eq.Q gives a relation between Ri, R^ and 7 of the ring when two Majorana states 
exist simultaneously, see the red line in fig.2 and Table I for 7 = 1/30 and Ri = 0.1. 

Table 1. External radius Re of the Majorana Ring with parameters Ri = 0.1^ and two 
values of 7. The radii are identified by the integer n. 

n 12 3 4 

7 = 1/30 0.256 0.531 0.742 0.9526 
7 = 1/15 0.559 0.987 1.411 1.835 

The superfluid density 

p{r) = r\Uir)f (10) 

of the two Majorana states for 7 = 1/30 and Ri = 0.1, Re = 0.953 are presented in Fig. 
3. One observes that due to exponentials in Eq.Q f± is located mostly near the internal 
(external) surface - the red and green lines respectively. 

Using the approximate periodicity of the Bessel functions in Eq.(|8]) with period of n one 
obtains for small 7 << 1: 

Re- Ri = 2TTn-f^ = Tm/2k± (11) 

where n is an integer. A finite system that conforms to these conditions will be denoted as 
"Majorana ring" in what follows. 



If the condition Eq.([9]) is violated (hence the rings will be termed "non- Major ana") 
the would be Majorana / = 1 states acquire a nonzero energy that oscillates around zero 
as function of Rf. for fixed Ri and 7. This is exemplified in Fig. 4 where energies of the 
/ = 1 states calculated numerically are given for 7 = 1/30 and Ri = 0.1 and a range of 
0.76^ < Re < 1.2^. The calculation utilizes the NAG Fortran Library Routine Document 
F02EBF. It computes all the eigenvalues, and optionally all the eigenvectors, of a real general 
matrix. One observes that at the "Majorana geometry" the energy of the electron branch 
(red line) changes sign. Similarly the hole branch has the opposite sign end the same absolute 
value of energy. This exhibits a phenomenon of " level crossing" in the BdG equation. Note 
that for a small sample and Re — Ri = vr (n + 1/2) /2A;_|_the energy of the / = 1 states might 
even exceed the gap magnitude A, see green lines in Fig. 4. 



IV. EXCITED STATES IN MAJORANA RINGS. 



For a vortex in an infinite s-wave superconductor, when the vortices are unpinned (freely 
moving) the low-lying spectrum of quasiparticle and hole excitations is equidistant, Ei = lu, 
where angular momentum / takes half-integer values. The "minigap" Emg in the s-wave 
superfiuids is of order oi u = A'^/Ep << A, where A is the energy gap and Ep is the 
Fermi energy. In the bulk chiral p-wave the spectrum of the low energy excitations remains 
equidistant, £'/ = (/ — 1) u, but now / is integer [15 1. 

For a Majorana ring there are excited states well below the continuum at / 7^ 1 typically 
close to the surface. The origin of their small energy is the vanishing of the superconducting 
gap on the nodes of p-wave superconductor. Such a state with minimal energy defines the 
"minigap" that protects the Majorana states from interferences due to the excitations. The 
excitation energies below the threshold (namely Andreev states) were calculated numerically 
for 7 = 1/30, Ri = 0.1 and Re = 0.98 , see Table 2. 



Table 2. Energy of Excited States in Majorana Ring with parameters Ri = 0.1^, i?e = 
0.953^,7 = 1/30. Energies of both the quasiparticles (-E+) and holes {E_) are given in units 
of A. 



6 



I 1 2 3 4 5 6 

E+ 2.6297-10-^ 0.1552 - 0.0135 0.5903 
E_ -2.6297- 10-^ - -0.0437 - - -0.1967 

Only the lowest (highest) energy E'^ (-£'-) for quasiparticles (holes) for each angular 
momentum is given (typically other excitations are beyond the threshold). At certain I 
there are no Andreev states. In the case of the small ring considered here one observes that 

the minigap Emg = 0.0135 appears at / = 4. 

The energy is just a fraction of A, still an order of magnitude larger than A'^/Ep for 
Sr2Ru04. The superfluid density p defined in Eq. flTU]) corresponding to its wave function is 
presented in Fig. 3 (blue line). The wave function even for such a small ring is peaked near the 
outer surface. To understand a relatively small value of the minigap compared to the one that 



protects the core Majorana states in the infinite system considered [18|], E^g = 0.25A , we 
have simulated finite large rings with Re = 10, 20. Of course the order parameter is no longer 
constant over such sizes and we have used A (r) = Atanh(r/^) for the order parameter 
distribution. One observes that there are surface modes with energies E (20^) = O.OIA, 
E (10^) = O.OIA. Therefore the energy of these modes is almost independent of Rg. 



V. STABILITY WITH RESPECT TO CHANGE OF THE ORDER PARAMETER 
DISTRIBUTION. 

The quantized geometry of the Majorana ring is not expected to be robust against pertur- 
bations like disorder, shape change, magnetic field distribution. We conjecture that varying 
parameters of the system, like the order parameter that depends on local temperature etc., 
one can still tune the system into a Majorana ring. To support this conjecture we calculate 
in perturbation theory the correction to the Majorana solution of the simple model Eq.([8]) 
in which the order parameter was approximated by a constant. Now we assume that the 
order parameter, in addition to a homogeneous rotation of the phase, depends on location 
s (r) = s {r, cp) = 1 + ip (r, cp). 

The only component of the BdG operator in Eq.([T]) that is corrected is the off 
diagonahL = - A {s (r) e^'^ (zV,. - Vy) +i [(^V, - Vy) {s (r) e*^)]} 
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L [s (r, ^)] 



A 



2kF 



^2i^ (d_ ^ %_d\ Lds_ _ \ ^ds_ 

2kF 1 r dip J r \ dip J dr 



(12) 
(13) 



The energy correction to the / = 1 Majorana states given in Eqs.Q to leading order in ijj is 



Ef, = (MH [s] - H f,) 







L[s]-L [1] \ Ub 



(14) 



Vb 



where a,b = ± (for the internal or external surface Majorana states of the unperturbed 
system). Using 



the diagonal elements are 



Et+ = Elf =1 I rU (r) e"^^ (L [s] - L[l]) [/+ (r) e'^] + cc = 0. 



(15) 



(16) 



The off diagonal read 

Em-" = Etr = [ ^ (r, ^) r/-/+ (r) (17) 

The condition for stability of the Majoranas therefore is: J^^, s (r, ip) rf_f^ = j^^ rf_f^. 
For any distribution that obeys this restriction, like for example s {r,ip) = 1 + S2 (v?) with 

S2 {(p) = 0, there is no splitting of the Majorana states. In particular this indicates that 
small deviations from the "Majorana ring" condition can be tuned away by controlling for 
example the value of the order parameter. 



VI. DISCUSSION AND CONCLUSIONS. 



To summarize, the spectrum of excitations of a single vortex in a chiral superconducting 
ring with internal and external radii Ri, Re comparable with the coherence length was 
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calculated. There is a pair of precisely zero energy states for Re — Ri = 7m/2k± for any 
integer n for k±^ > 1 where k± = sJ^m^Ep /fi^ — m±/miik'^ is the momentum in direction 
perpendicular to the magnetic field. Therefore a certain combination of geometrical factors, 
order parameter (dependent on temperature and material parameters) and magnetic field 
distribution make a "Majorana ring". They are not protected by topology, but are stable 
under certain deformations of the system. The condition has a character of a resonance 
phenomenon. 

The pair of exact Majorana states constitutes a qubit and might offer a promising 
method of the fault-tolerant quantum computation Uj. It is separated by a minigap 
of order A/70 from charged Andreev surface states, still order of magnitude larger than 
the Caroli,Matricon, deGennes minigap js] A'^/Ep for the superconductor Sr2Ru04 with 
A = 2K, Ep = lO^K. These states are generally ignored in other qubit proposals based 
on pair of Majorana states that belong to two different vortices. This configuration was 
comprehensively studied by Matsishita and Machida under assumption that the dis- 
tance between the vortices L is much larger than ^, so that tunneling between them can be 
treated as a perturbation. They found, using variational state made of the core Majorana 
states of each of the vortices, that the splitting energy of the two nearly Majorana states 
oscillates E± ~ Ai/fcpLcos (kpL + vr/4) e~^/^. Therefore at L = vr (n + 1/4) /kp one would 
get a degeneracy. This is reminiscent of our formula with two important differences. First 
our splittings are generally of order A, not Ay/kpLe~^^^ and second, more importantly, the 
Majorana states highly overlap, so that n is small, while for two vortices n > Lkp » 1. 

Practical proposals will depend on the thickness of the film. In the 2D limit k± = kp 
(the condition is determined below) and the Majorana condition is Re — Ri = Tin/2kp. For 
Sr2RuOi with j^, ^ = Q5nm the ring width should be a quantized in it /2kp = Anm ^ .^/30 
with minigap of order of IK. Moreover deviations from the Majorana condition lead to a 
very sharp splitting of order A, see Fig. 3. When the thickness of the film is D new channels 
open for the Andreev bound states: fc,, = 2?™,,//}, where n,, is integer. The new condition 
on the width becomes: Re — Ri = ^ [l — a (27m„/ Dk*)'^'^ ^ , where a = m±/m,, = 0.03 
is anisotropy and k* = \/2m_}^Ep/h. The second channel n^^ = 1 will enter for Re — Ri ^ ^ 
at width D = 27ra^/^/A;* = 5nm. Therefore the film is 2D only when it is thinner than 
that, while for thicker films, when several new channels are open, one of them can harbor 
Majorana states. 
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Figure captions 



Fig.l 

A p-wave superconducting mesoscopic disk with internal and external radii Ri and Re 
subject to a magnetic field perpendicular to the disk. 
Fig.2 

The resonance condition for the Majorana disk. The function x fo^' two values of 
the parameter 7 = 1/2^A;_l, 7 = 1/30 (blue), and 7 = 1/15 (red). To obey the resonance 
condition the internal and external radii Ri and Re should obey x {Ri) = X {Re)- 

Fig.3 

Superfiuid density (defined in Eq. fllOp ) of low energy states in Majorana ring with internal 
and external radii Ri = 0.1^ and Re = 0.953^ respectively and 7 = 1/30. The two Majorana 
/ = 1 states are the internal or core state (red) and external surface state (green) and the 
lowest energy excited charged state with / = 4 (blue). The excited state clearly resembles 
the surface states of p-wave superconductors. 

Fig.4. 

Violation of the resonance condition. The Majorana states are split into positive and neg- 
ative charged branches. The splitting energy in mesoscopic samples (like the one presented 
with Ri = 0.1^ and Re in the range 0.7^ — 1.3^) can be as large as the superconducting 
gap (green line) at Re = 0.8. One observes that for geometries far from the Majorana ring 
condition, Re = 0.742, 0.953 (for 7 = 1/30). Note the oscillation of splitting energy of 
quasiparticle and hole. 
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